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DESIGN  OF  NEAR-OPTIMAL  GUIDANCE  LAW 
FOR  SPACE  INTERCEPTION 

Shi  Xiaoping 

Simulation  Center, 

Harbin  Industrial  University 

ABSTRACT:  The  article  applies  an  augmented-linearization 
method  for  nonlinear  systems  to  propose  a  near-optimal  control 
resolution  method  of  secondary  type  index  problems  for  nonlinear 
systems.  This  method  is  applied  to  design  a  near-optimal 
terminal  guidance  law  for  space  interception.  As  indicated  in 
the  simulation  results,  the  near-optimal  guidance  has  good 
performance . 

Key  Words:  1.  terminal  guidance,  2.  optimal  guidance  law. 

I .  Introduction 

Classical  optimal  control  theory  has  been  relatively 
perfect;  however,  many  difficulties  will  be  confronted  when  using 
these  theories  to  solve  some  concrete  optimal  control  problems . 

In  particular,  when  a  system  is  nonlinear,  almost  no  analytical 
form  of  optimal  control  can  be  attained.  Therefore,  it  is 
necessary  to  seek  the  applicable  optimal  control  solution  method 
for  nonlinear  systems . 

Based  on  the  foregoing  reasons,  the  article  studies  the 
optimal  control  problems  of  a  type  simulation  nonlinear  system  in 
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the  secondary  index  significance.  By  applying  the  augmented- 
linearization  method  for  nonlinear  systems,  a  near-optimal 
control  solution  method  is  proposed.  As  an  application  example, 
this  near-optimal  control  solution  was  applied  to  design  a 
terminal-guidance  law  in  space  interception.  With  simulation 
comparison  against  the  optimal  guidance  law,  the  results  show 
good  effects. 

II.  Suggestions  of  the  Problem 

Assume  a  nonlinear  control  system 

.  x(t)=f(X(t))+BU(t),X(to)=*X. 

(2—1) 

In  the  equation,  X6R“'.U  eR“‘.T  6Ct.,t,]cR,B  ,  B  is  an  nxm- 

dimensional  constant  matrix;  f  is  a  continuous  function  of  x. 
There  exists  an  N+1  order  continuous  partial  derivative  with 
respect  to  Xi  (i=l,2, . . . ,n)  for  f.  N  is  a  fixed  natural  number, 
and  it  is  assumed  that  f(0)=0.  For  these  systems  of  f(0)  not 
equal  to  0,  it  is  very  easy  to  transform  to  f(0)=0. 

Now  it  is  required  to  find  the  optimal  control  vector  U(t) 
so  that  the  index 

J  =  jX^  (t,)SX  (t,)  +  yLcx^  (t  )Q  X  (t )  + 

U’'(t)RU(t)]dt  (2—2) 

is  the  smallest.  S,  Q,  and  R  are,  respectively,  the  steady 
weighted  matrixes  of  nxn  dimension,  nxn  dimension,  and  mxm 
dimension.  Moreover,  S  and  Q  are  nonnegative  steady  matrices, 
and  R  is  a  positive  steady  matrix. 

III.  Augmented  Linearization  of  Nonlinear  Systems 

By  introducing  the  row  vectors  of  1  to  N  order  monomials  not 
related  to  all  linearization,  as  derived  from  vector 


2 


X=[XiX2.  .  .X^]^  from  elements  Xi  (i=l,  2,  n), 

q4Cx,  ,X,,-  ,x,*,X,  ,x,,  -  ,X,'' ,  ••• , 
X/l^GR*^'  (3-1) 


N 

=  V  /-’I 


In  the  equation,  L=2:cI,+,_,,  and  assume  U=-KQ  (3-2) 


In  the  equation,  IkGR'“‘‘  /  then  Eq.  (2-1)  becomes 


X=f(X)-BKQ^g(K  ,X) 


(3-3) 


By  deriving  both  sides  of  Eq.  (3-1)  with  respect  to  t,  we  obtain 

q  =  [x,,X2-.nx.'"-’x„7 
[|gi  (K  »X)  .g2(K  *X)  1  •••  > 

NX.''-*g.(K  ,X)7 

4G(K,X)  (3-4) 


Since  the  Taylor  formula  applies  Gfc(K,X)  (k=l,  2,  ...,  L),  there 


G.(K.X)=G,(K.O)+xl[x.-|-f- 
1-1  jj 'ax, 

^yo,(K.o) 

■^(nTT), 


5j^7‘"‘G,(K,eX) 


(3-5) 


In  the  equation,  8  6(0. 1)CR 
written  as 


,  therefore  Eq.  (3-4)  can  be 


in  the  equation. 


Q(t)  =  b(K)+.i^(K)O(t)-60.N+,(O  (3-6) 


0N+.(l)  — 

V  1  ,  g.-^--*-'.G(K  ,X) 

.N.ri  r, !  — rj  *  aX'l-aX’S 

a 

x'-jnx''(t) 
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In  the  equation,  t  =  ex  •  Because  when  X=0,  £2=0,  and  f(0)=0, 
therefore,  g(K,0)=f  (0)-BK'0=0.  Then  from  Eq.  (3-4)  we  can  see 
Gk(K,0)=0,  (k=l,2,...,  L) .  Therefore,  from  Eq.  (3-5)  we  know 
that  b(K)=0  in  Eq.  (3-6).  Therefore,  we  have 

Q(t)=A(K)Q(t)  +  0N+i(t> 

If  we  neglect  the  term  /  then 

Z(t)  =  A(K)Z(t),Z(to)  =  £2(to) 

(3-8) 

In  the  equation,  zCR*"*'  •  The  foregoing  equation  is  the 

augmented  linear  system  of  system  (2-1)  obtained  with  N-order 
Taylor  expansion. 

IV.  Augmented  Linearized  Error  Analysis 

There  is  a  relationship  between  the  status  vectors  of 
nonlinear  system  (2-1)  and  augmented  linear  system  (3-8) 

X(t)  =  [l.„IO„a-,.]Z(t)  (4-1) 

Subtract  Eq.  (3-8)  from  Eq.  (3-7),  and  integrate,  then  we  obtain 

£3(t)-Z(t)  =  l'..H(t,T)0N+,(x)dT  (4-2) 

In  the  equation,  the  status  transfer  matrix  H(t,T)  satisfies 

|^^^i^=A(K)H(t,T)  (4-3a) 

iH(t,t)=I  (4-3b) 

Therefore,  we  can  obtain  the  status  approximate  error 

|X,(t)-Z,(t)KD(t)  liXir^’  (4-4) 

In  the  equation,  D(t)  is 


4 


>=  r  <ffH  ( t,  T )  r 
Jl.  ^  r,+.-+r.-N+l 


1  Ji+“+'.G(k  ,50 

r, !  -r. :  '  ax;,'*-2K„'" 


However,  in  the  equation  ||  x  11  is  defined  as 


li  =Sup|^ 


|Xi(t)I..  'i 


<Fr  =  C0-O,'O-OD 


Assume  that  D(t)  has  a  boundary  in  the  convex  domains  including  X 
and  0,  then  we  have 

lx,(t)-z,(t)l=o(  11  X  11  •‘'+00  =  1, 


2,—  ,n) 


(4-5) 


V.  Near-optimal  Control  of  Nonlinear  Systems 


Eq.  (3-8)  can  always  be  written  as: 

Z(t)=A(K)Z(t)-BKZ(t)  (5-1) 

In  the  equation,  [pj  •  Since  Z  is  an  approximation  of  fi, 
Eq.  (3-2)  can  be  approximated  as 


U  =  -KZ 


(5-2) 


Substitute  the  foregoing  equation  in  Eq.  (5-1)  and  we  obtain 

Z(t)=A(K)Z(t)+BU(t)  (5-3) 

In  the  equation,  zeR^'  ,  and  2(to)  =  [X,  (to),  — .X? 

now  the  performance  index  (2-2)  can  be  made  equivalent  to 

J  =  yZ^t,)SZ(t,)4-|-J  ..[Z^'t)QZ(t) 

-f-U^(t)RU(t)]dt  (5-4) 

In  the  equation 
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Therefore,  the  near-optimal  control  of  the  original- nonlinear 
equation  is 

u  =  — R"‘B^PZ  (5-7) 


VI.  In  designing  a  terminal-guidance  law  for  space  interception, 
the  space  interception  dynamics  equation  in  the  longitudinal- 
direction  plane  is  [1] 


r=V 


(6- la) 


V=ru5 
2V&)  a 


(6-lb) 

(6-lc) 


In  the  equations,  r  indicates  the  relative  distance  between 
interceptor  and  target;  V  stands  for  the  relative  velocity;  w 
indicates  the  angular  rate  of  the  line-of -sight;  and  a  indicates 
the  dynamic  acceleration  of  the  interceptor.  From  the  foregoing 
equation,  with  some  manipulations,  we  obtain 

(r  V  wO  ■- (  3  W + 2  V  coa  ) 

(6-2) 

From  this  equation,  by  introducing  the  status  variable  and  the 

control  variable,  X(t)=rVco*  (6— 3a) 

u(t)=3W-r2V&)a  (6-3b) 
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Then,  Eq.  (6-2)  becomes  a  first-order  nonlinear  system. 


In  the  equation,  V  can  be  considered  as  a  constant  because  in  the 
terminal  guidance  phase  of  space  interception,  generally  the 
orbital  parameters  of  target  motion  are  known.  The  interceptor 
aligns  with  the  target  in  flying  along  the  predeteinnined 
trajectory,  and  the  relative  velocity  between  interceptor  and 
target  varies  little  during  the  entire  terminal -guidance  phase; 
fundamentally,  this  is  a  known  constant. 

The  basis  for  designing  the  terminal -guidance  law  is  as 
follows:  let  the  angular  rate  e 

of  the  line  of  sight  approach  zero,  and  gradually  approach  the 
guidance  in  parallel.  In  this  article,  the  guidance  law  is 
designed  in  this  manner:  first,  an  index  is  presented  on  the 
system  (6-4) 

J=YSxHt,)+yL 

CqX*(t)+^u*(t)]dt  (6—5) 

Then,  the  optimal  control  u(t)  is  derived  so  that  J - >min.  Then 

from  Eq.  (6-3b),  the  guidance  law  is  derived.  Thus,  we  have 


In  the  following,  we  solve  for  the  problems  (6-4)  and  (6-5) 
for  near-optimal  control.  Here,  we  take  N=2.  In  other  words, 
the  system  (6-4)  is  augmented  linearized  as 


z= 


ID 


k  4--^^ 
2k  vj 


(6-6) 


In  the  equation,  Z=[x,x^]^.  Since 
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Z=[x,x=7.i? 

u(t)  =  — KZ=— k.x  — lijX^  (6-7) 


Therefore,  Eq.  (6-6)  can  be  written  as 

Z=A(K  )Z-f-Bu 

In  the  equation,  the  coefficient  matrix 


(5-8) 


A  (K  )  = 


O 

V2 

LD  2kJ 


Lo  J’ 


By  solving  for  the  Riccati  equation  [2], 

■-P(t)=P(t)A(K)+A^(K)P(t)- 
.P(t)B— i’‘P(t)+Q 


In  the  equation. 


P(t,)=s 


r<i  oi  [5  01 

Lo  oj  ’  Lo  oj 


(6-9) 


the  derived  P(t)  is  substituted  into 


K  =  [-—  0]P(t) 

(2 


Finally,  we  obtain 


--If 


VqjrSh(^ 

^t)+S-Ch(^ 

it) 

^'^ch(  J 

£t)+s  •sh(,&) 

(6- 10a) 
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kj  — 


2k?p-q 


VHq-4k?n) 

VT-ki 

2Vn-N/T-2ki  VIT) 

v^+k,-y7r 
■  2V*(vT+2k,>/r 

(6-lOb) 


icp( 


)t 


Substitute  the  two  foregoing  equations  in  Eg.  (6-7)  and  obtain 
the  near-optimal  control  of  the  problems  (6-4)  and  (6-5). 
Further,  based  on  Eg.  (6-3b)  we  can  derive  the  near-optimal 
guidance  law 

a  (t )  =  - 1- V to  (t )  -  (t )  CO  (t ) - 


!^r(t)vto*(t) 


(6-11) 


VII.  Nximerical  Simulation 

By  using  the  following  parameters  as  an  example  for 
numerical  simulations  S=0.4,  q=0.4,  n=0.2,  r(0)=150km, 

V=-7.5km/s,  ©(0 )=5xl0'^rad/s .  Then  the  initial  condition  of 
system  (6-4)  is  X(0)=-281.25(m^/s^  [sic]).  Figs.  1  and  2 
indicate  the  solution  curves  of  optimal  control  problems 
(6-4)  and  (6-5).  Solid  lines  represent  the  near-optimal  control 
derived  by  using  the  method  in  the  article;  the  dashed  lines 
indicate  the  optimal  control  derived  from  the  numerical  method  by 
using  the  NPSOL  software  package.  Curves  in  the  figures  indicate 
that  the  near-optimal  control  is  relatively  close  to  the  optimal 
control . 

Figs .  3  to  6  indicate  the  comparative  relationship  between 
the  near-optimal  guidance  law  (represented  by  solid  lines) 
derived  by  using  the  method  in  this  article,  and  the 
corresponding  optimal  guidance  law  (represented  by  the  dashed 
lines)  derived  by  using  the  NPSOL  software  package. 
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From  these  figures,  we  can  see  that  they  are  relatively 
close  to  each  between  the  optimal  guidance  law  and  the  near- 
optimal  guidance  law.  Additionally,  two  physical  quantities  of 
relative  distance  and  relative  velocity  of  interceptor  and  target 
exhibit  basically  no  difference  under  optimal  control  and  near- 
optimal  control . 

The  article  was  received  for  publication  on  August  4,  1994. 


Fig .  1 .  Optimal  and  near-optimal  status 


Fig.  2.  Optimal  and  near-optimal  control 
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Fig.  3.  Relative  distance  between  interceptor  and  target 


-31 

^  I 

3 

i 

>  -si 


— '-1 

I 

-Ti 


Fig.  4.  Relative  velocity  between  interceptor  and  target 
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Fig.  5.  Variation  of  angular  rate  of  line  of  sight 
Rate  of  LOS  indicates  the  angular  rate  of  line  of 
sight 
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Fig.  6.  Variation  of  powered  acceleration 
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